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ABSTRACT 

We consider a gamma-ray burst (GRB) model based on an anisotropic fireball with an axisymmetric 
energy distribution of the form e(9) oc 9~ k , and allow for the observer’s viewing direction being at an 
arbitrary angle 9 V with respect to the jet axis. This model can reproduce the key features expected from 
the conventional on-axis uniform jet models, with the novelty that the achromatic break time in the 
broadband afterglow lightcurves corresponds to the epoch when the relativistic beaming angle is equal 
to the viewing angle, 9 V , rather than to the jet half opening angle, 9j. If all the GRB fireballs have such 
a similar energy distribution form with 1.5 < k< 2, GRBs may be modeled by a quasi-universal beaming 
configuration, and an approximately standard energy reservoir. The conclusion also holds for some other 
forms of angular energy distributions, such as the Gaussian function. 

Subject headings: gamma rays: bursts - shock waves - ISM: jets and outflows 


1. INTRODUCTION 

Recently, several independent approaches have led to 
the conclusion that long gamma-ray bursts (GRBs) have 
a standard energy reservoir of several 10 50 ergs (Frail et al. 
2001, hereafter F01; Panaitescu & Kumar 2001, hereafter 
PK01; Piran et al. 2001, hereafter P01). An important 
ingredient of this argument is that the putative jet open¬ 
ing angles, 9j, as inferred from the afterglow lightcurve 
breaking times, tb, have a broad distribution, but just of 
the right form to compensate for the wide dispersion of 
the “isotropic” energy emitted in y-rays, E Jt i so , so that 
E 1 = (E 7i i so /47r)(0|/2) is essentially invariant (F01). The 
total energy of the fireball should be E to t > E 7 + Eq, 
where Eq is the initial kinetic energy of the fireball in 
the afterglow phase assuming an adiabatic evolution, and 
the inequality takes into account the possible energy loss 
during the radiative regime in the early afterglow phase 
that has evaded the present observations, as well as en¬ 
ergy losses outside the y-ray band (e.g. the BATSE win¬ 
dow) or in non-electromagnetic forms (e.g. neutrinos and 
gravitational waves) during the prompt phase. Writing 
E 1 = r]Et 0 t where ry is the gamma-ray emission efficiency, 
Utot could be mainly contributed by Eq if rj is small (e.g. 
< 0.1). PK01 and P01 found that Eq is also distributed 
in a narrow range. For a uniform jet, this leads to the 
inference that Eo = (dE/dCl)(9 2 /2) is also essentially in¬ 
variant. However, in the above analysis, and in the current 
afterglow jet models which are used to determine 9j , it is 
generally assumed that the jets are uniform, with sharp 
cut-offs at the edges, and that the line-of-sight cuts right 
across the jet axis. None of these assumptions are neces¬ 
sarily true in general (Meszaros , Rees & Wijers 1998; Mac- 
Fadyen & Woosley 1999; Woods & Loeb 1999; Nakamura 
1999; Paczyriski 2001; Salmonson 2001; Dai & Gou 2001). 
On the other hand, although it is not difficult to construct 
a central engine model which makes GRBs with a standard 
energy reservoir but with quite different beaming angles, 
it would be more elegant to have a model that all the GRB 
beams share a standard energy reservoir as well as a quasi- 


universal beaming configuration (M. J. Rees, 2001, private 
communication). Here we show that such a model can be 
constructed by taking account of the off-axis anisotropic 
jet effects, without violating the present observational con¬ 
straints. 

2. THE MODEL 

Our assumption is that all the long GRBs have a quasi- 
universal beam configuration, with a strong anisotropy of 
the angular distribution of the fireball energy around an 
axial symmetry. The jet axis is physically related to the 
rotational axis of the central engine, so it is reasonable to 
assume that initially the closer to the jet axis, the higher 
the energy concentration. The actual angular distribution 
of the fireball energy is unknown, and we model it as (e.g. 
Meszaros et al. 1998) 

dE/dfl = e(0, </>) = e(0) = e 0 9~ k , (1) 

within the range 9 m < 9 < 0, where 9,„ is a very small 
angle within which some deviation from (|]J) is necessary to 
avoid the divergence at 9 = 0, and 0 is some large angle 
which exceeds the presently measured 0 3 by at least a fac¬ 
tor of two (for the simplification of the discussions below). 
The real angular energy distribution may differ from the 
power law (jlj), but most of our discussions below can be 
generalized to other forms of distribution functions (e.g. 
see Q and relevant discussions below). The adoption of 
(jT) is for the simplicity of the discussions. The angular- 
dependence of the baryon loading rate is uncertain, and 
we assume that it is weak so that the Lorentz factor angu¬ 
lar distribution follows a similar law, i.e., r(0) oc 9~ k (of 
course, the law should be modified when r(0) approaches 
unity). We make furthermore the assumption that 

(1.5) < k <2 (2) 

in ( 0 )- The reason for this requirement will become evi¬ 
dent later. The main conjecture of the model is that the 
dispersion in the afterglow data of the breaking time, tb, 
is a manifestation of the diversity of viewing angles of the 
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observers, rather than to the diversity of intrinsic opening 
angles of the jets themselves. In other words, that what 
were inferred by Frail et al. (2001) as 0j are essentially 
0 V in our model, where 0 V is the observer’s viewing angle 
with respect to the jet axis. We will test whether the above 
hypothesis is able to pass the following three criteria: (i) 
When r(0„) 1/0 V , the jet dynamics along the line-of- 

sight satisfies the isotropic law T(9 v ,t) oc t~ 3 / 8 (for sim¬ 
plicity, we only discuss an adiabatic fireball running into 
an interstellar medium with a constant density), where t 
is the observer time, and r(0„,f) is an effective Lorentz 
factor assuming an isotropic fireball which could mimic 
the emission in the direction 6 V at the time t ; (ii) When 
r(6g<l/0„, the dynamics changes so that the lightcurves 
steepen; (iii) The total jet energy E to t is essentially a uni¬ 
versal value. 

For an isotropic adiabatic fireball running into a uniform 
medium, T(f) oc e 1 ^n~ 1 ^{at)~ 3 ‘^ , where e = dE/dfl is 
the energy per solid angle, n is the ambient medium num¬ 
ber density, and the blastwave radius is written in a gen¬ 
eral form as R = aT 2 ct, where the factor a effectively takes 
into account the surface of equal-arrival-time as well as the 
thickness of the emitting region. To test the criteria (i), 
the key is to estimate the effective energy per solid angle, 
e(9 v ,4> v ,t ) = e{9 v ,t), in the direction (9 v ,<f> v ), and to eval¬ 
uate the possible time-dependence of this value. When 
r(^,t)=r>l/#„, the observer can only observe a solid 
angle around (9 v ,tf v ) with a half opening angle of order 
1/F due to the relativistic beaming effect 1 . By definition, 
the effective energy per solid angle in the direction ( 9 V , <f> v ) 
is 


e(9 v , (f v , i) — e(9 v , t) 


lUv-i/r s i n 0d9 


rOv + i/r 
Je v - l/r 


sin 9d9 


( 3 ) 


due to the axial symmetry. In the small angle approxima¬ 
tion, which is relevant to the present discussions 2 , one has 
sin# ~ 9. When F 1/9 V and noticing ( 5 ]), this gives 


Ie V+ l/r e o 01 kd 0 

e(0v,t) * v ~/ +1/r -^ to9f k = e(9 v ). 

ItBd0 


( 4 ) 


This is a time-independent quantity, since the sideways 
expansion effect is not important at the same stage (see 
discussions below). We then have 

f(#„,i) oc [e(#„, t)] 1 / 8 n _1 / 8 (at) -3 / 8 oc t -3 / 8 , T 1 /9 V . 

( 5 ) 

This indicates that the observer does not feel the 
anisotropy of the fireball when the relativistic beaming 
angle 1/T is much smaller than the viewing angle 9 Vl but 
observes the fireball as if it were isotropic. This is the same 
conclusion as drawn in the on-axis uniform jet model. The 
conclusion (|lj) does not require and holds for any k 
value. In fact, it even holds for some other forms of angu¬ 
lar energy distributions, e.g. the Gaussian distribution, 

e{9) = e 0 exp[—(l/2)(6»/6» 0 ) 2 ], (6) 


as long as the first order Taylor expansion term of these 
functions lead to the same result. Notice that the factor a 
of various forms may deviate from the conventional value 
(e.g. ~ 4), since the shape of the equal-arrival-time sur¬ 
face will be distorted due to the anisotropic distribution of 
the fireball energy. However, its time-dependence, if any, 
would be very small. Therefore, although it may influ¬ 
ence the absolute values of the afterglow flux levels, such 
an effect does not change the blastwave dynamics in the 
viewing direction. 

In principle, the jet configuration is time-dependent, due 
to the effects such as the energy redistribution and the lat¬ 
eral expansion. In the lab frame, the causally connected re¬ 
gion subtends an angle of A 9 ~ c s t'/R = c s /cT{9 ), where 
t' = R/cT{9) is the comoving time since the explosion 
along the 9 direction, c s is the expansion speed, which may 
be either the relativistic speed of sound ~ c/y/ 3 (Rhoads 
1999), or simply the speed of light (Sari et al. 1999). As 
r(ftu) 1/0 V , the line-of-sight direction is causally discon¬ 
nected from other regions, so the dynamics evolves essen¬ 
tially independently, so that the description in (|3|) holds. 
When the blastwave decelerates so that the line-of-sight 
bulk Lorentz factor r(#„) drops close to and below 1/9 V , 
the dynamics along the line-of-sight starts to change, and 
r(0„) will deviate from the oc t -3 / 8 dependence. There 
are several effects that play a role. First, as 1/T(9 V ) ex¬ 
ceeds 9 V , the observer starts to feel the energy deficit due 
to the drop of the energy distribution, i.e. deviation of the 
power law i), on the other side of the jet axis. Although 
the calculation of e(9 v , t) is no longer straightforward, this 
deficit effect should mimic that in the uniform jet model 
as long as k is not too flat, say, k > 1.5. Second, the 
anisotropic jet has a trend to resume the isotropic shape. 
As r(0„) ~ 1 /9 V , the viewing direction starts to connect 
the jet axis causally. The energy outflow from the cone 
defined by 9 V becomes prominent, and this equivalently 
decreases e(9 v ) in the viewing direction. In the meantime, 
the initial material within the 9 V cone starts to spread 
into a wider cone, and the observer would feel a stronger 
deceleration, although the global sideways expansion will 
become evident only when F(#„) drops below 1/0. All 
these effects tend to steepen the afterglow lightcurve, al¬ 
though the degree of steepening is unclear without detailed 
numerical dynamical calculations. In any case, in the 
asymptotic phase of sideway expansions, R is essentially a 
constant, and one would eventually have T(9 v ,t) oc f -1 / 2 
(from t ~ R/ r 2 , Rhoads 1997, 1999; Sari et al. 1999). 
In this regime, the temporal indices of the lightcurves in 
the various spectral regimes would follow closely the same 
predictions as in the uniform jet models (Sari et al. 1999; 
Rhoads 1999). For example, in the slow-cooling regime 
(which is usually the case after the viewing or “jet” break), 
for spectral regimes both below and above the cooling fre¬ 
quency, the asymptotic spectral flux is F v oc t~ p , where p 
is the power-law index of the electron number distribution. 
For reasonable values of p (e.g. ~ 2.2), this is consistent 
with several GRB afterglow observations. The above dis¬ 
cussion should also hold for other distribution functions 


1 Strictly speaking, the observer will see a smaller half cone on the close side of the jet axis, and a larger half cone on the far side to the 

8 's, due to different Lorentz factors in different directions. This will modify the integral limits in (H). but does not influence the conclusion in 
and the relevant discussions. 

2 The largest “jet” angle in F01 is 0.411, and the approximation is good within 3%. 
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such as (|]), mainly because eventually all the initial con¬ 
figurations will be smeared out. However, to address the 
lightcurves properly within different models, including the 
relevant gradual transition between asymptotic regimes, a 
detailed dynamical description and numerical calculation 
is necessary, and we postpone this to a future work. 

We have shown that the present model can reproduce 
the key features of the on-axis uniform jet model, with an 
arbitrary k value as long as it is not too flat. The next 
question is whether the model can also retain the merit 
of a standard energy reservoir invoked in the conventional 
jet model. In principle, one does not have to fulfill this 
constraint, but just wishes so for the sake of elegance. By 
definition, the total energy in a fireball with an energy 
distribution given by (Jl]) is 

/* © /*© 

E tot = 2tt / e(0) sin 9d9 ~ 27r / eg9 1 ~ k d9. (7) 

Jo Je m 

For k < 2 and 0>0 m , we get 

0-7r 

E tot ^ g —b 2 ~ k e(O v )e 2 v , (8) 

where we have parameterized 0 = bd v . We can see that 
the quantity 2ne(9 v )6% (which is essentially the E 1 of F01, 
or Eg of PK01 and P01) is quasi-invariant, if k and E to t 
are constant (or have a small scatter). The only extra 
scatter is introduced through the scatter of b, which is in¬ 
troduced by the scatter of 9 V (assuming the same 0 for 
all GRBs). However, for the index (2 — k) this scatter is 
greatly reduced if k is not much smaller than 2. This is 
another reason why we require, say, k > 1.5, in (|2j). A 
smaller 0 can also reduce the b scatter. Notice that the 
b scatter tends to raise E to t in GRBs with smaller 0 v ’s 
(and hence larger b’s), which seems to be helpful to re¬ 
duce the Eg scatter in PK01. An important implication 
of equation (||) satisfying such a constraint is that the to¬ 
tal energy reservoir is standard, but the absolute value 
need no longer necessarily be several times 10 50 ergs, but 
would depend on the value of k and the typical value of b. 
Given reasonable values, .Etot could be one order of mag¬ 
nitude higher than that of F01 and PK01, but this could 
be still well accommodated within conventional central en¬ 
gine models (Meszaros , Rees & Wijers 1999). The closer 
k approaches 2, the larger the standard energy reservoir 
one requires. At k = 2, equation (||) should be modified 
in a form containing a logarithmic term, and the energy 
requirement is the highest (see discussions in Rossi, Laz- 
zati & Rees 2001). Also the scatter of 6 m must be very 
small for k = 2, while for k < 2, the actual value of 9 m is 
not important. For k > 2, generally E to t (eq.Q) can not 
be expressed in terms of e(9 v )9%, since most ofthe energy 
is distributed at small angles. The standard energy bud¬ 
get argument no longer holds 3 . A quasi-universal beaming 
configuration as well as a standard energy reservoir is how¬ 
ever in general obtained if the requirement (]2;) is satisfied 
for an energy distribution such as ©• 

For other forms of energy distributions, a standard en¬ 
ergy reservoir is also attainable. For example, for the 
Gaussian distribution (|6|), one has E to t ~ eg&g. 


3. DISCUSSION 

We have shown that an off-axis anisotropic jet with an 
energy distribution with angle given by equation (|p (or 
other forms such as ||), is able to reproduce the key ob¬ 
servational features of a conventional on-axis uniform jet 
model, e.g. such as producing a “jet break” signature in 
the light curve. The novelty here is that the achromatic 
break time tb in the broadband afterglow lightcurves no 
longer corresponds to the time when the relativistic beam¬ 
ing angle is equal to the jet half opening angle, 9j. Rather, 
it corresponds to the time when the relativistic beaming 
angle is roughly equal to the observer’s viewing angle 9 V 
relative to the jet axis. In this model, the broad distribu¬ 
tion of tb in the data is no longer due to the intrinsic scatter 
of the jet opening angles among different bursts, but is at¬ 
tributed to the distribution of the observer’s lines of sight. 
For a power law energy distribution ([I]) (or a Gaussian en¬ 
ergy distribution |J), if the constraint ||) is satisfied, all 
the GRBs may have a quasi-universal beaming configura¬ 
tion, besides a quasi-standard energy reservoir. We deem 
this to be a more elegant picture than the conventional 
on-axis uniform jet model. In addition, the homogeneous 
nature of the conventional model is more idealized, and 
the present inhomogeneous model is likely to be a closer 
representation of what could be expected in nature. 

The predictions of this inhomogeneous model for the af¬ 
terglow lightcurves are not completely equivalent to those 
of the uniform jet model. The key difference should oc¬ 
cur around the “jet break” time. Our model should give 
a more gradual variation at the break than the uniform 
jet model, which assumes a sharp drop off at the jet edge. 
The so far sparsely studied sideways expansion effect in 
an anisotropic jet may further complicate the problem. 
The shape of the break should also depend on the angular 
energy distribution function and some unknown param¬ 
eters, such as k. Detailed modeling is necessary in or¬ 
der to address these questions. In any case, the gradual 
break expected in our model is not inconsistent with sev¬ 
eral well studied afterglow lightcurves, and some simula¬ 
tions have shown that the conventional jet models usually 
also give gradual and smooth jet breaks (e.g. Panaitescu 
& Meszaros , 1999; Moderski, Sikora & Bulik 2000; Huang 
et al. 2000). Both models are compatible with the present 
data, but this situation may change as better data becomes 
available and as more detailed simulations are performed. 

Recently, Rossi et al. (2001) have independently dis¬ 
cussed the power-law model (jlj) in more detail. They plot¬ 
ted the afterglow lightcurves for the k = 2 case which 
mimic those of the on-axis uniform jet model, and also 
discussed the more general cases of k ^ 2. Here we have 
presented a general analytical argument, showing that the 
blastwave dynamics at the line of sight is identical to the 
uniform jet model in the asymptotic regime for a locus of 
models of the general form of equation (pj) , as long as k is 
not too flat. With this particular form of the angular de¬ 
pendence of the energy, in order to have a standard, finite 
energy reservoir for all bursts one requires the constraint 
(|2). The upper end k< 2 of the constraint (|2]) ensures that 
the total energy can be expressed in terms of e(9 v )9% and 


3 However, if 0 m is not too small (e.g. a not very small fraction of 9 V ), the case A:>2 could still retain the feature of a standard, finite (but 
even larger) energy reservoir. 
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does not diverge. (However, the case k >2 could also have 
the same virtue if 9 m is not too small compared with 9 V ). 
The lower end of the constraint, k > 1.5, ensures that 
the scatter introduced by 9 V is not too large, and that the 
energy-deficit effect at the other side of the jet is not too 
small. We have also found that the main features in the 
power-law model are also applicable to some other angular 
energy distributions, e.g. such as the Gaussian form ([|). 

For ease of discussion, we have here assumed that the 
upper limit of validity of the assumed angular distribution 
is 0 > 2 9 V . This is to avoid that the observer feels the 
energy deficit beyond 0 before the relativistic beaming 
angle exceeds 9 V . Indeed, if 0 < 2 0 V , T(9 V ) starts to devi¬ 
ate from the value predicted by the adiabatic law oc f -3 / 8 
after it is less than (0 — 9 V )~ 1 . In this regime, the upper 
limits for ^-integration in the numerators of both (||j) and 
(Q) should be replaced by 0. Thus the maximum correc¬ 
tion factor with respect to the 0 > 2 9 V case is a factor of 
0/2 9 V . Even for 0 = 9 V (i.e., the line of sight marginally 
cuts the jet edge), the deviation is at most a factor of 1/2. 
We therefore conclude that the 0 effect may in most cases 
not be important. The main reason is that the large angles 
contribute a small portion of the total energy in the beam 
due to the distribution of the form (jlj). 

In our model, the “isotropic” luminosity function will be 
determined by the assumed angular distribution, N(e)de = 
N{9)d9 oc sin 9d9 oc 9d9 (the latter being for small 9). 
From equations ( 0 ), (§ and substituting e by L, we get 
the luminosity function predictions in our model, e.g., 

N(L)dL oc L~ 1 ~ 2 ^ k dL (9) 

for the power-law model, and 

N(L)dL oc L~ 1 dL (10) 

for the Gaussian model. To test these luminosity func¬ 
tions, redshift measurements are needed. Using only the 
bursts for which optical redshifts have been determined so 
far, e.g. as compiled in F01, PK01, the above luminosity 


distributions are not consistent. However, this discrep¬ 
ancy could be due to small number statistics (~ 20 in all 
or ~ 10 on each side of the mid-point). The fact that the 
small sample size or other selection effects related to the 
afterglow detections could lead to a spurious inconsistency 
is also suggested, for example, by the clear deficit of low lu¬ 
minosity (e.g. possibly due to large viewing angle) bursts 
at higher redshifts (z > 1) in FOl’s data set. Alternatively, 
some other distance indicators have been proposed, such 
as spectral time-lags, e.g. Norris, Marani & Bonnell 2000, 
or variability measures, e.g. Fenimore & Ramirez-Ruiz 
2002, Reichart et al. 2001 (and interpretations of these 
indicators in terms of the viewing angle have been dis¬ 
cussed by, e.g., Salmonson & Galama 2002 and Norris et 
al 2002). If one accepts such distance indicators and their 
inferred redshifts at face value, the observational GRB lu¬ 
minosity function inferred for a much larger bursts sample 
(e.g. Schaefer, Deng & Band 2001) is not inconsistent with 
the theoretical distribution |)- And, using a different ap¬ 
proach, Schmidt (2001) obtained a flatter observational 
GRB luminosity function, which over a large range of lu¬ 
minosities is compatible with the model distribution (fLOj). 

A natural consequence of this model is that the dis¬ 
tribution of break times tb, and hence the 9 V distribu¬ 
tion, should be related to the statistical distribution of 
viewing angles and to the shape of the beam distribution. 
The present data and the preliminary calculations are not 
sufficient to draw firm constraints on parameters of such 
models. However, the comparison of such predictions or 
more detailed versions of them against future data in the 
Swift era should provide interesting constraints, as a larger 
quantity of more accurate redshift measurements become 
available. 
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ulated this research, to E. Rossi for sending us their paper, 
and to B. Paczynski, J. Granot, E. Waxman, L. J. Gou for 
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by NASA (NAG5-9192 and NAG5-9153). 
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